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Abstract 

Wc investigate the optimal configurations of n points on the unit sphere for a 
class of potential functions. In particular, we characterize these optimal config- 
urations in terms of their approximation properties within frame theory. Fur- 
thermore, we consider similar optimal configurations in terms of random distri- 
butions of points on the sphere. In this probabilistic setting, we characterize 
these optimal distributions by means of special classes of probabilistic frames. 
Our work also indicates some connections between statistical shape analysis and 
frame theory. 
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1. Introduction 

Frames are overcomplete (or redundant) sets of vectors that serve to faith- 
fully represent signals. They were introduced in 1952 by Dufhn and Schacffer 
12 1, and reemerged with the advent of wavelets 0, [H, 15, 18, 21 1. Though the 



overcompleteness of frames precludes signals from having unique representation 
in the frame expansions, it is, in fact, the driving force behind the use of frames 
in signal processing d, [2^ [2^ . 

In the finite dimensional setting, frames are exactly spanning sets. However, 
many applications require "custom-built" frames that possess additional prop- 
erties which are dictated by these applications. As a result, the construction 
of frames with prescribed structures has been actively pursued. For instance, 
a special class called finite unit norm tight frames (FUNTFs) that provide a 
Parseval-type representation very similar to orthonormal bases, has been cus- 
tomized to model data transmissions |y,bO| . Since then the characterization and 
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construction of FUNTFs and some of their generalizations have received a lot of 
attention djUMEll- Beyond their use in applications, FUNTFs are also related 
to some deep open problems in pure mathematics such as the Kadison-Singer 
conjecture [vf. FUNTFs appear also in statistics where, for instance, Tyler 
used them to construct M-estimators of multivariate scatter [s^] ■ We elaborate 
more on the connection between the M-estimators and FUNTFs in Remark 12.31 
These M-estimators were subsequently used to construct maximum likelihood 
estimators for the the wrapped Cauchy distribution on the circle in 2J| and for 
the angular central Gaussian distribution on the sphere in (ssj . 

FUNTFs are exactly the minimizers of a functional called the frame poten- 
tial [H. This was extended to characterize all finite tight frames in [s^. Fur- 
thermore, in 1^ 22 1, finite tight frames with a convolutional structure, which 



can be used to model filter banks, have been characterized as minimizers of an 
appropriate potential. All these potentials are connected to other functional 
whose extremals have long been investigated in various settings. We refer to 
[13, for details and related results. 

In the present paper, we study objects beyond both FUNTFs and the frame 
potential. In fact, we consider a family of functionals, the p-frame potentials, 
which are defined on sets {xi}^i of unit vectors in Mf^; see Section[51 These po- 
tentials have been studied in the context of spherical i-designs for even integers 
p, cf. Seidel in and their minimizers are not just FUNTFs but FUNTFs 
that inherit additional properties and structure. Common FUNTFs are recov- 
ered only forn = 2. In the process, we extend Seidel's results on spherical 
t-designs in [29| to the entire range of positive real p. 

In Section [3l we give lower estimates on the p-frame potentials, and prove 
that in certain cases their minimizers are FUNTFs, which possess additional 
properties and structure. In particular, if < p < 2, we completely characterize 
the minimizers of the p-frame potentials when N = kd for some positive integer 
k. Moreover, when iV = d-l- 1 and < p < 2, we characterize the minimizers of 
the p-frame potentials, under a technical condition, which, we have only been 
able to establish when d = 2. We conjecture that this technical condition holds 
when d > 2. Finally in Section 21 we introduce probabilistic p-frames that 
generalize the concepts of frames and p-frames. We characterize the minimizers 
of probabilistic p-frame potentials in terms of probabilistic p-frames. The latter 
problem is solved completely for < p < 2, and for all even integers p. In 
particular, these last results generalize [1^ as well as the recently introduced 
notion of the probabilistic frame potential in [l6l |. 

Fiirther relations to statistics: Besides the results on FUNTFs used in 
[iSjISls^, and mentioned above, frame theory has essentially evolved indepen- 
dently of statistical fields such as statistical shape analysis [IJ] and directional 
statistics Nevertheless, there still exist several overlaps, and to the best of 
our knowledge, these overlaps have not yet been fully explored. Recently, frame 



theory has been used in directional statistics [17[, where FUNTFs are utilized 
to investigate on statistical tests for directional uniformity and to model and 
analyze patterns found in granular rod experiments. We must point out that 
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similar results were obtained earlier by Tyler in [33j . 

Probabilistic tight frames are multivariate probability distributions whose 
second moments' matrix is a multiple of the identity, and they are used in [lij to 
obtain approximate FUNTFs. The latter approximation procedure is connected 
to a classical problem in multivariate statistics, namely estimating the popu- 
lation covariance from a sample, which is closely related to the M-estimators 
addressed in [24 , 32 , 3^ . The p- frame potentials and their probabilistic counter- 
parts that we consider in the sequel, are linked to the notion of shape measure, 
shape space, and mean shape used in statistical shape analysis. In Section [521 
we establish a precise connection between the full Procrustes estimate of mean 
shape [2^ Definition 3.3] which is the eigenvector corresponding to the largest 
eigenvalue of the frame operator. Moreover, this eigenvalue coincides with the 
upper frame redundancy as introduced in [j]. The full Procrustes estimate of 
mean shape also saturates the upper frame inequality. Moreover, the p-frame 
potentials form size measures as required in statistical shape analysis, and their 
minimizers among all collections of N points on the sphere define a shape space 
modulo rotations. 

We hope that the present paper will renew interests in more investigation 
on the role of frames and the p-th frame potential in directional statistics and 
statistical shape analysis. 



2. The p-frame potential 

2.1. Background on frames and the frame potential 

To introduce frames and their elementary properties, wc follow the textbook 

i- 

Definition 2.1. A collection of points {xi}^^ C M'^ is called a finite frame for 
R'' if there are two constants Q < A < B such that 

N 

A\\xf <^\{x,x,)f < B\\x\\^, for all X e R''. (1) 

i=l 

If the frame bounds A and B are equal, the frame {xi}^i C M'' is called a finite 
tight frame for M.'^. In this case, 

N 

=^\{x,x^)f, foran.TGM''. (2) 

i=l 

A finite tight frame {xi}fL^ C M'' consisting of unit norm vectors is called a 
finite unit norm tight frame (FUNTF) for . In this case, the frame bound is 
A = N/d. 

A collection of unit vectors {xi]^^^ C S'^~^ is called equiangular if there 
exists a nonnegative constant C such that \{xi, Xj) \ = C, for i ^ j. 
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Given a collection of N points {xi}^i in W^, the analysis operator is the 
mapping 

Its adjoint operator is called the synthesis operator and given by 

JV 

Using these operators, it is easy to see that {xi}fLi is a frame if and only if the 
frame operator defined by 

N 

S -.^ F*F -.R"^ ^R'^, x^"^{x,x^)xi 

1=1 

is positive, self-adjoint, and invcrtiblc. In this case, the following reconstruction 
formula holds 

N 

X = ^(S'"^Xi,x)x,, for all x e M"*, (3) 

and {S~^Xi}fLi, in fact, is a frame too, called the canonical dual frame. If 
{xi}fLi is a frame, then {S~^^^Xi}fLi is a finite tight frame. Moreover, note 
that {xi}^i is a FUNTF if and only if its frame operator 5 is ^ times the 
identity. 

As mentioned in the introduction, the question of the existence and char- 
acterization of FUNTFs was settled in 0] , where the frame potential, defined 

by 

N N 

FP({xJiIi)=5]5]|(x„x,)p, (4) 

i=l 3 = 1 

was introduced and used to give a characterization of its minimizers in terms of 
FUNTFs. More specifically, they prove the following result: 

Theorem 2.2. 0, Theorem 7.1] Let N be fixed and consider the minimization 
of the frame potential among all collections of N points on the sphere S'^~^ . 

a) If N < d, then the minimum of the frame potential is N . The minimizers 
are exactly the orthonormal systems for R'' with N elements. 

b) If N > d, then the minimum of the frame potential is . The minimizers 
are exactly the FUNTFs for M'' with N elements. 

We shall prove in the sequel that the frame potential is just an example in 
a family of functionals defined on points on the sphere, and whose minimizers 
have approximation properties similar to those of the frame potential. But 
first, we briefly comment on the relation between FUNTFs and M-estimators 
of multivariate scatter: 
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Remark 2.3. The concept of FUNTFs is used in signal processing to represent 
a signal a; G R'' by means of a; = X]i!=i (^5 ^i)^*' which is similar to the 
well-known expansion in an orthonormal basis. FUNTFs have also been used in 
statistics to derive M-estimators of multivariate scatter: For a sample {xi}^^ C 
Mf^, where Xi 0, for i = 1 . . . , N, Tyler aims to find a symmetric positive 
definite matrix F such that 

M(T) = - V 

^tt ^^^^^ 

is the identity matrix. Whenever this is possible, the estimate V of the popula- 
tion scatter matrix is then given by V ~ T^^ . We refer to 3^ for details. Note 
that M (F) = Id implies that 

i- ' Xt \ I i ' a;,, I ^ ^d-i 



i=l 



forms a FUNTF. Moreover, {x,}fL^ C S'^~^ is a FUNTF if and only if M{Id) = 
Id. 

2.2. Definition of the p- frame potential 

Definition 2.4. Let be a positive integer, and < p < oo. Given a collection 
of unit vectors {xi}^^ C S"^^^, the p-frame potential is the functional 

N 

When, p = oo, the definition reduces to 

FPoo,Af({a;i}iIi) = sup |(a;,;,.Tj)|. 

It is clear that the p-frame potential generalizes the frame potential in (|4]). 
Finite frames and the p-frame potential also extend to complex {zi}^^ C 
S^^^ — {z £ : \\z\\ = 1} and are related to statistical shape analysis, a 
tool to quantitatively track the physical deformation of objects. We refer to 



14 Chapters 2, 3 & 4] for more details on shape analysis, but we briefly indi- 
cate here its link to the p-frame potential. The shape of an object is specifled 
by landmark points that altogether form the shape space. Often, a suitable 
transformation is applied flrst in order to study shape independently on the 
object's size. To remove the feature of size, we must specify a size measure g 
( \l4 . Definition 2.2]), which is a positive function defined on S^~^ that satisfies 

= sg{{z,}fL,), for all s G R+ and {z,}fL, C S^'- 
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It is immediate that the following family of functionals can be seen as size 
measures: 

In particular, gi represents the centroid size (when the shape is centered at 
the origin), which is one of the most common size measures in statistical shape 
analysis. Given two complex configurations Zi, Z2 € C'^ derived from landmarks 
that code two-dimensional shape, the full Procrustes distance ([3, Definition 
3.2]) is 

-I3e - a-ib 



dFizi,Z2) = inf 

Given TV configurations {zi}fLi C 
shape is defined by 



l^il 
ird 



Z2\ 



the full Procrustes estimate of mean 



N 



zp = arg \j\i_^{^d],{z,,z)), 



cf. IJ, Definition 3.3]. The average axis from the complex Bingham maximum 



likelihood estimator is the same as the full Procrustes estimate of mean shape 
for two-dimensional shapes, and the same holds for the complex Watson dis- 
tribution, cf. [13, Sections 6.2 & 6.3]. If we assume that {zi\^^-^ are centered 
around zero, then zp is given by the eigenvector corresponding to the largest 
eigenvalue A of the frame operator of the normalized collection {-jif^lfli. Fur- 
thermore, one observes that this eigenvalue A is the upper frame redundancy of 
{zi\f^i as introduced in Q, which also coincides with the optimal upper frame 
bound B in (HI). Therefore, the full Procrustes estimate of mean shape satisfies 
the upper frame inequality with equality. Moreover, we have observed that the 
root mean square of the full Procrustes estimate of mean shape is 1 — . 

Before stating the next elementary result, we recall some basic definitions in 
physics: A conservative force F, is a vector field defined on R'', such that —F is 
the gradient of some potential P that is then induced by the conservative force. 
Lemma 12.51 below was proved for the frame potential in and we extend it to 
all p- frame potentials with 1 < p < 00: 

Lemma 2.5. For each p e (l,oo), the p-frame potential FPpjv is induced by 
the conservative force Fp = fp{\\o- — 6||)(a — b), for a,b G S"^"^ , where 

\p{l - for \x\ < V2, 

l^Pi^ — 1)^^^, otherwise. 



fp{x) := 



Fp is a central force between the 'particles' a and b that we call the p-frame 
force. 



Proof. The function 



Pp[x) 




for \x\ < V2, 
otherwise, 
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is differentiable and satisfies Pp{x) = —xfp{x). This is sufficient to verify 
that the potential Pp{a,b) := Pp{\\a — b\\), defined for a,b & 5"*"^, satisfies 
WaPpici, b) = —F{a, b), where b is held fixed. Thus, Fp is a conservative vector 
field. The physical meaningful potential Pp{a, b) is in fact given by 

Pp(a,6) = Pp(||a- fell) = |l - i||a - 6f |^ = |(a,6)|^ 

where we used that ||a|| = ||6|| = 1. Consequently, the p-frame potential is 
induced by the conservative central force Fp. □ 

As a consequence of the above lemma, {xi}^^^ C S'^~^ are in equilibrium 
under the p-franic force if they minimize the p-frame potential among all collec- 
tions of points on the sphere. Note that such a collection of equilibria modulo 
rotations form a shape space. 

Remark 2.6. Wc will use repeatedly the fact that for a fixed {x^j^^ C S''^"^, 
the p-framc potential ¥¥p^N{{xi}f^i) is a decreasing and continuous function 
of p e (0,03). 



3. Lower estimates for the p-frame potential 

We start this section with a few elementary results about the minimizers of 
the p-frame potential as well as their connection to t-designs. In fact, potentials 
on the sphere, and t-designs have been well investigated [H [13, [3, [H, 34 1. 
However, one of the key differences between t-designs and our p-frame potential 
is that the former is considered only for positive integers t while the latter is 
investigated for p S (0, oo). 

3.1. The Welch bound revisited 



If p = 2fc is an even integer, one can use Welch's results [36| to conclude 
that, for {xj^i C S"*-!, 

FPp.^({x.}f=i) > ppy. (6) 

We shall verify that this estimate is not optimal for small A^, by proving an 
estimate for FPp jv when 2 < p < oo. The following Proposition first appeared 
in [ill: 

Proposition 3.1. Let {xj^^ C N > d, and 2 < p < oo, then 

FPp,^({xJf^i) > NiN^l){^^Y^' + N, (7) 
and equality holds if and only if {xi}fLi is an equiangular FUNTF. 
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Proof. For 5 = ^ + Holder's inequality yields 

U{^^,^,)WJ\\i. < \\i{^u^,))^^MN{N-l)y/^. (8) 

Raising to the p-th power and applying 7 = 5 ~ ^ leads to 

< \\i{x.,X,)).,^,\\l{N{N-l)r/'-\ (9) 

Therefore, 

Using the fact that I (^^i, a^j) P > -^^ N (see Theorem 12. 2p implies that 

whieh proves ([7]). 

To establish the last part of the Proposition, we recall that an equiangular 
FUNTF {xk}k=i C satisfies 



= y for ah i 7^ j (10) 

see, [HI, Isij, for details. Consequently, if {xk}^^i is an equiangular FUNTF, 
then ([7]) holds with equality. 

On the other hand, if equality holds in ([7]), then X^j^j ~ ^ ^ 

and {xijfLi is a FUNTF due to Theorem 12.21 Moreover, the Holder estimate 
([8|) must have been an equality which means that a;j)| — C for i ^ j, and 
some constant C > 0. Thus, the FUNTF must be equiangular. □ 

By comparing ([6]) with ([7]), it is easily seen that the Welch bound is not 
optimal for small N: 

Proposition 3.2. Let {xj^^ C S"'"^ and p = 2k > 2 be an even integer. If 
d<N < C+J;"^), then 

FPp,^({^4f=i) > m - +N> (11) 

Proof. The condition on N implies 1 > ^d+k-i'^ , and adding (A^ — 1) ( ^(jvl'^i) ) > 
to the right hand side leads to 

Multiplication by N and Proposition 13.11 then yield (jlip . □ 
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Remark 3.3. The estimate in Proposition l3.1l is sharp if and only if an equian- 
gular FUNTF exists. In 3l|, Sections 4 & 6], construction (and hence existence) 
of equiangular FUNTFs was established when d + 2 < < 100. For general 
d and N, a necessary condition for existence of equiangular FUNTFs is given, 
and it is conjectured that the conditions are sufficient as well. The authors 
essentially provide on upper bound on N that depends on the dimension d. 
Therefore, Proposition 13.11 might not be optimal when the redundancy N/d is 
much larger than 1. 

3.2. Relations to spherical t-designs 

A spherical t-design is a finite subset {xi]^^i of the unit sphere S'^~^ in M'', 
such that, 

N 



^ V/i(a;.0 = / h{x)da{x), 



for all homogeneous polynomials h of total degree equals or less than i in d 
variables and where a denotes the uniform surface measure on S'^~^ normalized 
to have mass one. The following result is due to 34, Theorem 8.1] (see [l^] 
for similar results). 



Theorem 3.4. fSA . Theorem 8.1] Let p ~ 2k be an even integer and {xi} 
{-x,}l, c 5^-1, th 



N 
i=l 



en 



T?v \N ^^ 1 ■ 3 ■ 5 ■ ■ ■ (p - 1) 2 

FP.Ml-.h^i) > d{d + 2)---{d + p-2f ' 

and equality holds if and only if {xi}^^i is a spherical p-design. 

3. 3. Optimal configurations for the p- frame potential 

We first use Theorem 12.21 to characterize the minimizcrs of the p- frame po- 
tential for < J) < 2 provided that the number of points is a multiple of the 
dimension d: 

Theorem 3.5. Let < p < 2 and assume that N = kd for some positive integer 
k. Then the minimizers of the p- frame potential are exactly the k copies of any 
orthonormal basis modulo multiplications by ±1. The minimum of ([S]), over all 
sets of N = kd unit norm vectors, is k^d. 

Proof. If we fix a collection of vectors {xi}fLi, then the frame potential is a 
decreasing function in p G (0,2). Therefore, 

N^d = FP2M{^^}ll) < FPpM{^^}ll)■ 

Now suppose that {yi}fLi consists of k copies of an orthonormal basis 
{ejf^ifor R'^. Then, 

minFPp,^({xJiIi) < FPp,^({yJiIi) = FP{{y,}l,) = N^d. 
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Consequently, 



minFPp.jv({a;af=i) = min FP({xfc}Li) = N^/d. 

Thus, FPp jv has the same mmimuni as PP. Clearly k copies of an orthonormal 
basis of W'- form a FUNTF, and hence minimize the FP due to Theorem 12.21 
Consequently, the k copies of an orthonormal basis also minimize the p-frame 
potential for < p < 2. On the other hand, the inner products must be or 1 
to obtain a minimizer. The smallest p-frame potential then have the k copies of 
an orthonormal basis. □ 

Wc now consider the p- frame potential for = + 1 . The case p = 2 is 
covered by Theorem 12.21 Note that any FUNTF with N = d + 1 vectors is 
equiangular 0,1331 • Hence, the case 2 < p < oo is settled by Proposition 13. H 
so we focus on p <E (0, 2). 

One easily verifies that, for po = \og(d) — ^' orthonormal basis plus one 
repeated vector and an equiangular FUNTF have the same po-frame potential 
FPpp^£i+i. Under the assumption that those two systems are exactly the mini- 
mizers of YVp^^^d+ii the next result will give a complete characterization of the 
minimizers of FPp^d+i, for < p < 2. However, we have only been able to 
establish the validity of this assumption when d = 2, cf. Corollarv 13.71 

Theorem 3.6. Let N ^ d + \ and set po = iog{d) = log(jv-i) ■ 
{xi}iLi C 5'''""'^, and assume that FPpg^N{{xi}^i) > N + 2, with equality holds 
if and only if {xi}^i is an orthonormal basis plus one repeated vector or an 
equiangular FUNTF. Then, 

(1) for < p < pq, then for any {xi}fLi C S'^^'^ , we have FPp^pf{{xi}^^) > 
N + 2, and equality holds if and only if {xi}^i is an orthonormal basis 
plus one repeated vector, 

(2) for Po < p < 2, then for any {xi}fLi C S'^~^, we have FPp^jv({a;,;}^i) > 
2po (Nd)^ PO +N, and equality holds if and only if {xi}fLi is an equian- 
gular FUNTF. 

Proof. Under the assumptions of the Theorem, let < p < pa, then 

minFPp„,Ar({.T,}iIi) < minFPp,Ar({a;0.=i), 

and using {yi}^i consisting of an orthonormal basis of M'^ with one repeated 
vector, yields that 

minFPp,jv({2:JiIi) < FPp,jv({2/,:}iIi) = iV + 2 = minFPp„,w({xj£i). 

Consequently, minFPpojv({a;i}^;^) = min FPp^jv({a;,;}i^i) ~ N + 2. Since an 
orthonormal basis plus one repeated vector minimizes the p-frame potential for 
p = Po, it must also minimize FPp.Ar for < p < po, which proves (1). 
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Assume now that po < p < 2. Choose r such that ^ = ^ + ^, Holder's 
inequality yields 

< \\{{x.,T,)),^,\\l{N'~Nr/^. 
Since we assume that (2) holds, we have ^ 2, which leads to 

> 2po {Nd) Po . 

This concludes the proof of (2). By applying ([10]), one then checks that an 
equiangular FUNTF satisfies (2) with equality. 

The "only" part comes from the fact that the Holder inequality becomes 
an equality only if the sequences are linearly dependent. This means that the 
{xi}iLi ^-re equiangular. They must then satisfy \{xi,Xj)\ = i. Thus, by (jlOp 
they form an equiangular FUNTF 

When d = 2 we can in fact verify the main hypothesis of Theorem [321 which 
leads to the following result: 

Corollary 3.7. Let {.xjf^i C S'^ , and set pa = jSSgl. Then, 

FPpo,3({xOLi) > 5, 

and equality holds if and only if {xi}f^^ is an orthonormal basis plus one repeated 
vector or an equiangular FUNTF. 
Consequently, 

(1) for < p < Pq, then for any {xi}'^^^ C , we have FPp^{{xi}^^i) — 
and equality holds if and only if {xi}^^i is an orthonormal basis plus one 
repeated vector, 

(2) for Po < p < oo, then for any {a;i}f^i C , we have FFp^^dxi}^^^) > 
^ +3, and equality holds if and only if {xi}f^^ is an equiangular FUNTF. 

The minimum of FPp.3, for < p < 00 is plotted in Figure [Tj 

Proof. Clearly (1) and (2) follows from Theorem 13.61 once the minimizers of 
FPpo,Af are characterized. 

Without loss of generality, let /? be the smallest angle between xi, X2, and 
X3 and let a be the second smallest angle between them. This yields, of course, 
0< (3 <a. 

Case 1 : For 0<a + /3< f,we have 

i E = cosP(a) + cosP(/3) + cosP(a + /3) = F(a,/3). 
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Figure 1: Minimum of FPp,3 in Corollary |3^ for < p < 10. 

Since I < po, the p- frame potential is differentiable in a and (3, and its critical 
points are 

= —pcos^~^ (a) sin(a) — pcos^~"'^(a + /3) sin(a + /?) 
= -pcosP-\(3)sm{f3) - pcosP~i(a + /3) sin(a + 

This imphes that either a=:/3 = Oor/3 = and a — ^. In the first case, we 
have a maximum since it implies xi = X2 = X3. The latter case means that 
two points are identical and the third one is perpendicular which is a potential 
minimum of the p-frame potential. 

Case 2: We can assume that j < a (otherwise we are in Case 1). We can 
further assume that f<a<|<a + ;3<^ {ii ^ < a + 13 < n. Otherwise, 
substitute Xi with —Xi. Wc now have 

i ^ I (x„ Xj)\P = coaPia) + cosP(/3) + (- cos(a + P))p = G(a, /3). 

The critical points of G are given by 

= -pcosP"^(a) sin(a) +p(-cos(a + I3))p^'^ sm{a + P) 
= -pcosP-^{l3) sin(/3) + p{~ cos(a + sin(a + /3). 

By subtracting one equation from the other and raising to the second power, 
we obtain 

z^"^(l - zi) = - ^2), 

where Zi = cos(a)^ and Z2 = cos^(a + (3). Since sin^(x) — 1/2 > cos^(a;) 
for all 7r/2 < x < 2tt/3, we have < zi < 1/2 and < Z2 < 1/4 because 

Z2<l-Z2- 1/2. 

We consider the function F{z) = zP~^{1 — z) on < z < 1/2. F achieves 
its maximum at z = w 0.3691 and is convex, cf. Figure [2J Therefore, 
for < zi < 1/2 and < Z2 < 1/4, we have F{zi) = F(z2) if and only if 
zi = Z2 or zi = 1/2 and Z2 = 1/4. For zi = cos^(q!) = 1/2, we have a = tt/S 
and Z2 = cos^(7r/3 + /3) = 1/4 yields f3 = n/3. The case zi = Z2 leads to 
cos^(q;) = cos^(q! + /3) which implies a + /3 — 7r/2 = tt/2 — a. This is equivalent 
to /? = TT — 2a. Since we assume /3 < a, we obtain tt/S < a < tt/2. 



12 



(a) = - 1), < ^ < I (b) /(a) = 2cosP(a) + {-cos{2a))P, f < a < 

Figure 2: F and / in the proof of Theorem 13.71 

We now check the minima of the function 

/(a) = cosP(a) + cosP(7r - 2a) + (- cos(7r - a)Y 
= 2cosf(a) + (-cos(2a))f, 

on 7r/3 < a < 7r/2, cf. Figured] Its derivative 
df 

-^(a) = -2pcosP-\a) sin(a) + 2p(- cos(2a))P~^ sin(2a) 
oa 

vanishes if and only if 

cos^^^(a) sin(a) = (sin^(a) — cos^(a))^~^2 cos(a) sin(Q;) 
For a 7^ 7r/2, this yields 

cosf-i(a) = (1 - 2cos2(a))f^i2cos(a). 
By substituting x = cos(a), we obtain, for < .t < 1/2 

xP-^ = (1 - 2x^)P-^2x, 

which is equivalent to 

1/2 = - 2x)P-^x <^0^{l/x- 2a;)a;i/(f-i) - (1/2)1/(p-^) 
and define the new function 

= -2x'+i - (l/2)«, 

where q = l/{p—l). To show that g has only one extremal point on < a: < 
we differentiate 

||(x)--2(q + l)x^ + (Q-l)x«-^ 
The term §f (a;) vanishes if and only if 

(g-l)x-«-2 =2(g + l)x'', 
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which is equivalent to 

2 ^ g-1 ^ 2^1og2(3) 
2{q+l) 21og2(3) ■ 

Hence, a; sa 0.3618 and ^ does not have any other zeros on < a; < 1/2. This 
means that g has only one extremal point and can then only have two zeros on 
< a; < 1/2. The zero at a; = 1/2 corresponds to a minimum. This means that 
the zero of If at a; « 0.3618 is a maximum of g. Hence the other zero of g is 
between and ~ 0.3618. However, this other zero corresponds to a maximum 
of /. The minimum of / can thus be at a; = or x = 1/2. This implies a — 7r/3 
or a = 7r/2. It is easy to verify that a — 7r/3 would lead to 13 — 7r/3, and 
a = n/2 yields /3 = 0. Thus, the minimum of the p- frame potential corresponds 
to either an orthonormal basis plus one repeated element {a — 7r/2, /3 — 0) or 
an equiangular FUNTF {a ~ (3 ~ 7r/3). One easily checks that both situations 
lead to the same global minimum. 

□ 

In view of Theorem 13.61 and Corollarv l3.7[ we have the following conjecture: 
Conjecture 3.8. Let N = d+l and po = l^^fl) ■ Then 

FPp„M{^k}Li)>N + 2, 

and equality holds if and only if{xi}^^ is an orthonormal basis plus one repeated 
vector or an equiangular FUNTF. 

One can check that 1 < po < 2, for d > 1. According to Proposition 13. 11 the 
minimizers of the p-frame potential for 2 < p < oo are exactly the equiangular 
FUNTFs. Thus, our conjecture essentially addresses the range < p < 2. 

Remark 3.9. Although we are primarily interested in real unit norm vectors 
{xi}iLi C M"^, we should mention that Theorem 12.21 Theorem 13. 5( Equation 
([6]), and Propositions 13.11 and 13.21 still hold for complex vectors C 
that have unit norm. The constraints on N and d that allow for the existence 
of a complex FUNTF are slightly weaker than in the real case (sij . 



4. The probabilistic p-frame potential 

The present section is dedicated to introducing a probabilistic version of 
the previous section. We shall consider probability distributions on the sphere 
rather than finite point sets. Let A4{S'^~^ , B) denote the collection of probability 
distributions on the sphere with respect to the Borel sigma algebra B. 

We begin by introducing the probabilistic p-framc which generalizes the 
notion of probabilistic frames introduced in jlil ]. 
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Definition 4.1. For < p < oo, we call ^ G 7^(5'^ ^, B) a probabilistic p-frame 
for M'* if and only if there are constants A, B > such that 

A\\yr<[ \{x,y)\Pdtiix)<B\\y\\P.. Vy G R"^. (12) 

We call /i a ii^/if probabilistic p-frame if and only if we can choose A ~ B. 

Due to Cauchy- Schwartz, the upper bound B always exists. Consequently, 
in order to check that fi is a probabilistic p-frame one only needs to focus on 
the lower bound A. 

Since the uniform surface measure a on S'^~^ is invariant under orthogonal 
transformations, one can easily check that it constitutes a tight probabilistic 
p-frame, for any < p < oo. 

Given a probability measure /i G B), we call 

F -.R"^ ^ Lp{S'^-\n), x^{x,-) 
the analysis operator. It is trivially seen that 

ll-F'(y)llLp(s<'-i,M) = \\{x,y)\\L^{S'^-\,,) < \\y\\ 

for all < p < oo. The dual of F is called synthesis operator and is given by 
F* : Lg{S''-\^,) ^ R", f^f f{x)xd^i{x), 

where 1 ^ - + -, and I < p < oo. In fact, F* is well-defined and bounded 
operator on all Lr{S'^~^ , fj.) where 1 < r < oo. Indeed, for / e LriS"^^^ , fj.) we 
have 

\\F*{f)\\< f \f{x)\\\xW{x)^ f \f{x)\df,{x)<\\f\\L^. 

Given /i € A4{S'^^^ , B), the second moments matrix of /i is the d x d matrix 
defined by 

S = F*F =(^j x,Xjdfi{x)^ (13) 

with respect to the canonical basis for K''. As we show next, the second moments 
matrix S plays a key role in determining if /i is a probabilistic p-frame. We refer 
to 0, where a result was proved for similar discrete p- frames. 

Proposition 4.2. //I <p < oo, then fj. e M{S'^^^,B) is a probabilistic p-frame 
if and only if F* is onto. 

Proof. As mentioned earlier the upper bound in the probabilistic p-frame defi- 
nition always holds. So we only need to show the equivalence between the lower 
bound and the surjectivity of F*. 
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Assume that F* is surjective. Since F is injective, then S is invertible and 
for each y € M*^ we have 

llyjp = \{Sy,S-'yM = \{F*Fy,S-'yM = \{Fy, FS'^y) ^,^^,,\, 

which can be estimated as follows: 

\\yf < ll^^(2/)l|L.|l^^(^"'y)!lL.' < C\\Fiy)U4S-\y)\W. < C\\Fiy)\\MR''- 
Therefore, for all y 7^ £ M'^ we have 1/C||y|| < \\F{y)\\LP that is 

l/C\\yr< I \{x,y)\^dii{x). 

Thus ^ is a probabilistic p-frame. 

Now assume that /i is a probabilistic p-frame, but that F* is not surjective. 
Then, there exists z 7^ G M'* such that {z,F*f) = for all / e Lp,{S'^-'^ , ^i). 
Consequently 

(z,F*/)k. = {F{z)J)l,-,l^, = / f{x){x,z),d^,{x) = 

for all / e Lp'. A contradiction argument leads to {x,z) — for all x G S"^^^ 
which implies that z = 0. Thus F* is surjective. □ 

The second moments matrix can also be used to show that a probabilistic 
p-frame gives rise to a reconstruction formula that extends the finite frame 
expansion in ([3]). In addition, the next result generalizes the reconstruction 
formula for tight probabilistic frames obtained in [igI . Lemma 3.7]. 

Proposition 4.3. Let Q < p < 00 and assume that p, G Ai{S'^~^,B) is a 
probabilistic p-frame for R'^. Set // = /i o S". Then for each y G Mf^ , we have 

y^j Sz{z,y) dfliz) = z {S z , y) djl{z) . (14) 

Proof. The result follows by noticing that y ~ SS^^y = S^^Sy. □ 
The above result motivates the following definition: 

Definition 4.4. Let < p < 00. If /i e M[S'''^^ , B) is a probabilistic p-frame 
with frame operator 5, then fl ^ p o S e M{S^^{S''-^^), S^^B) is called the 
probabilistic canonical dual frame of fi 

Note that if n in Proposition 14.31 is the counting measure corresponding to 
a FUNTF {xi}fL^, then fi is the counting measure associated to the canonical 
dual frame of {a;i}^i. 

Lemma 4.5. a) If p is probabilistic frame, then it is a probabilistic p-frame for 
all 1 < p < 00. Conversely, if fj. is a probabilistic p-frame for some 1 < p < 00, 
then it is a probabilistic frame. 

b) Let 1 <p < 00. If fj. is a probabilistic p-frame, then so is the canonical dual 
/Et. 
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Proof, a) Assume that /i is a probabilistic frame and let 1 < p < oo. Then, we 
only need to check that the lower inequality of (|T2|) holds, since the correspond- 
ing upper bound is trivial. By Proposition 14.21 (applied to p — 2), S = F*F is 
invertible and for each y gR'^ wc have 

\\yf = \{Sy,S-'yM = \{F*Fy,S-'yM = \{Fy, FS-'y) ^,^^,,1 

which can be estimated as follows: 

hr < \\m\\LjFiS-'y)h^, < C\\Fiy)UJ\S-\y)U. < C\\F{y)hJy\U.. 

Therefore, for all y ^ e R'^ we have 1/C||y|| < \\F{y)\\L^ that is 

l/C\\yr< f \{x,y)\Pdt,{x). 

For the converse, assume that /i is a probabilistic p- frame for somep > 2. Then, 
for all y 7^ e M'^, 

A\\yr< f \{x,y)\Pdt,{x) 

^1 \{x,y)\^\{x,y)r^dii{x) 
< I \\xr-^\\yr-^\{x,y)\'d^,{x) 
= \\y\r' I \{x,y)\H^,{x), 
from which it follows that 

A\\yr< f \{x,y)\'dfi{x). 

If /i is a probabilistic p-framc for some p <2. Then, for all y 7^ e M'', 
\\yf = \{Sy,S-'y)^A = \{F*Fy,S-'yM = \{Fy,FS-^y) l^^lJ, 
which can be estimated by 

hf <\\Fy\\L^FS~'y\\L^, <C\\Fy\\LMl 

where we have used the fact that for p < 2, L2{S'^~^ , fi) C Lp{S'^~'^ , fi). This 
conclude the proof of a). 

b) If /i is a probabilistic p- frame for some 1 < p < 00, then by a) /i is a 
probabilistic frame. In this case, fl is known to be a probabilistic frame, cf. [l6j . 
and thus a probabilistic p-frame. □ 

We are particularly interested in tight probabilistic p- frame potentials, which 
we seek to characterize in terms of minimizers of appropriate potentials. This 
motivates the following definition: 
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Definition 4.6. For < p < oo and yu G M.{S'^ ^t^), the probabilistic p-fn 
potential is defined by 



PFP(/.,p)=/ / \{x,y)\Pd^lix)d^liy). (15) 

From tiie weak;-star-compactncss of tiie collection of all probability distribu- 
tions on the sphere, we can deduce that PFP{fi,p) admits a niinimizer which 
satisfies 

PFP(n) = min PFP(u,n). (16) 

We now turn to the minimizers of the probabilistic frame potential PFP(/z). 
In the process, we extend some ideas developed in Q to the probabilistic frame 
potential. 

Proposition 4.7. Let < p < oo and let /i be a minimizer of psp . then 

(1) Jg,_, \{x,y)\Pdtiix) = PFP(p), for all y G supp(m), 

(2) Jg,_, \{x,y)\Pdt,{x) > PFP(p), for all y G S''-\ 

Proof. The proof will use the following observation. Let /i be a probability 
measure on S'''"^ and choose a measure v, such that i'{S'^^'^) = and fi+eiy > 0, 
for all < e < 1. Let us also introduce the notation 

PFP(^,z.,p) := / / \{x,y)\Pd^lix)dl^iy). 

We then obtain 

PFP(^) < PFP{n + ei^,p) 

= PFP(/i,p) PFP(z^,p) + 2ePFP(/^, i^,p) 
= PFP(/i) PFP(i^,p) +2ePFP(^,i/,p). 

We thus have < ePFP(z^,p) + 2 PFP(/i, i/,p), for all < e < 1, which implies 
PFF{n,v,p) > 0. 

We now prove (1) using a contradiction argument. In particular, assume 
that (1) does not hold. This implies that there are yi,y2 G supp(/i) such that 



\{x,y2)\Pd^iix) < / \{x,y,)\Pd^iix) =: b. 

Set P^{y) = /gd-i \{^Ty)Y'dii{x). Let K be an open ball around yi in S''^^ and 
so small that y2 ^ K and that the oscillation of {y) on K is smaller than . 
Let m = > 0- One can check that the measure v defined by 

v{E) ■.= m5y^{E) - ^i{Er\K), E e B, 
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satisfies v{S'^ ^) = 0, and ^ + ev > Q. Hence, PFP(/^, z/,p) > 0. On the other 
hand, we can estimate 

VYF{^i,v,p)^ [ P^{v)du(y) = P^{y2)m^ j P^iy)dfi{y) = am- [ P^(2/)dAi(y) 



and so 



b — a b ~ a 
PFP(/i, v,p) < am — (o — )m = — m < 0. 



This is a contradiction to PFP(/i, v,p) > and implies that there is a con- 
stant C such that P^{y) = C, for all y e supp(/i). We still have to verify that 
the constant C is in fact PFP(p): 

PFP(p) = m{n,p) = f P^{y)dii{v) 

Lipp(/i) 

CdM(y) = C. 

The proof of (2) is similar to the one above, and so we omit it. □ 
The following result is an immediate consequence of Proposition 221 
Corollary 4.8. Let < p < oo and let ^ be a minimizer of p5|) . then 

(1) \{x,y)\Pdt,{x) = PFFip)\\y\\P, for all ^ G supp(/i), 

(2) supp(/i) is a complete subset o/M^. 

Proof. (1) directly follows from (1) in Proposition l4.7l 

(2) If supp(/i) is not complete in M'', then there is an element y G S'^^^ that 
is in the orthogonal complement. However, this contradicts (2) in Proposition 
K7\ □ 



We can now characterize the minimizers of the probabilistic p- frame potential 
when < p < 2. In fact, we shall show that these minimizers are discrete 
probability measures, and the following theorem is the analogue of Proposition 

KB 



Theorem 4.9. Let < p < 2, then the minimizers of (|15p are exactly those 
probability distributions p that satisfy both, 

(i) there is an orthonormal basis {xi, . . . , Xd] for such that 
{xi, . . . , Xd} C supp(/^) C {±xi, . . . , ±Xd), 
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(ii) there is f : S"^ ^ ^ M. such that fJ.{x) = f{x)i'±xi,...,±xd{^) 

The measure zyj-^j ^-^.^(.t) in Theorem 14.91 denotes the comiting measure of 
the set {zLxi : i = 1, . . . ,d}. 

Prool Since < |(x,y)| < 1, for x,y e S"^'^, we have PFP(/x,p) > PFP(/x,2). 
In [la . Theorem 3.10] it was shown that the normalized counting measure 
2^X1, ...,xa of f-ii orthonormal basis minimizes PFP(-, 2). Due to PFP(i^':ri,...,2:d, 2) 
VYY'{^Vxl....^xa^P)^ obtain that 2^x1,. ..,xa also minimizes PFP(-,p) and hence 
PFP(p) = PFP(2). 

In the following, we prove that all minimizers of PFP(-,p) are essentially 
induced by an orthonormal basis. Let /i be a minimizer and let v,w € supp(/i). 
We first show that w)| <E {0, 1}. The implications = 1 if and only if 

V = w and (u, w) = — 1 if and only if v = —w arc trivial. 

Suppose now that v ^ ±w and {v,w) 7^ 0, then there exist e > and > 
such that 

(a) Be{v) n Be{w) = and ^l{Beiv)), ^l{B,{w)) > S^. 

(b) for aU x e Bg{v) and y e B^{w), |(a;,y)|^ > \{x,y)\'^ + e. 
By using B = B^{v) x B^{w), this implies 

PFP(/x,p)= / \{x,y)\Pdfiix)dfi{y) + f \{x,y)\Pdfiix)dfi{y) 

JB J S''-^yS''-^\B 

> I {\{x,y)\^ +e)dfi{x)dfi{y)+ / \{x,y)\^dti{x)dfi{y) 

JB J S-^-^xS-^-^XB 

= PFP(m, 2) + e^iiB,iv))^iiB,iw)) 

> PFP(/^, 2) + £,52 > PFP(/x, 2), 

which is a contradiction. Thus, we have verified that |(a;, y)| G {0,1}, for all 
x,y G supp(/i). Distinct elements in supp(/x) are then either orthogonal to each 
other or antipodes. According to Corohary 14. 8[ supp(/i) is complete in R'^. 
Thus, there must be an orthonormal basis {xi}f^i such that 

{xi,...,Xd} C supp(/x) C {±a;i, . . . , ±2;^}. 

Consequently, there is a density / : S'^^^ — > R that vanishes on S'^^^ \ supp(//) 
such that ^(.t) = /(a;)i^±xi,...,±a;d(2;)- 

To verify that / satisfies (ii), let us define / : S"^^^ ^ R by 



fi^) + f{~x), X e {xi,. . . ,Xd} 

0, otherwise. 
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This implies that fl{x) = f{x)i'xi,...,xa{x) is also a minimizcr of PFP(-,2). But 
the minimizers of the probabihstic frame potential for p = 2 have been inves- 
tigated in [3 Section 3]. We can follow the arguments given there to obtain 
/(xi) = i, for aU i = 1, . . . ,d. □ 

For even integers p, we can give the minimum of PFP(/x,p) and characterize 
its minimizers. The following theorem generalizes Thcorcm l3.4l Moreover, note 
that the bounds are now sharp, i.e., for any even integer p, there is a probabilistic 
tight p-frame: 

Theorem 4.10. Let p be an even integer. For any probability distribution /i on 
PFP(m,p)=/ / \{x,y)\Pdti{x)dfi{y)> ^ ' ^ ' 5 • • • (p - 1) 



d-1 J qd-1 



d{d + 2)---{d+p-2)' 

and equality holds if and only if fi is a probabilistic tight p-frame. 

Proof. Let a = | — 1 and consider the Gegenbauer polynomials {C"}n>o defined 
by 

Co"(x) = l, C'^{x) = 2ax, 
CZ{x) = h^2x{n + a- l)C^i(x) -{n + 2a~ 2)C'^_^{x)] 

{C^'^Yn=i is an orthogonal basis for the collection of polynomials of degree less 
or equal to s on the interval [—1,1] with respect to the weight 

i.e., for m ^ n, 

CirHx)Ci:\x)wix)dx = 0. 

'1 



They are normalized by 



Ci"\x) il-x^)"-^dx 



n!(n + a)[r(a)]2 ' 

The polynomials t^, p an even integer, can be represented by means of 

p 



k=0 

It is known (see, e.g., [H [l^l) that > 0, i = 0, . . . ,p, and Aq is given by 



1 



c 



Ao = - / tPwit)dt, 



-1 
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where 

7r2'^+3r(d - 2) 



c = 



(|-i)r(|-i)2 

Moreover, induces a positive kernel, i.e., for {xi}'^^-^ C S"^^^ and {ui}^^-^ C 
R, 

JV 

^ >0, VA; = 0, 1, 2, ... 

ij = l 

see [l|, [3] ■ Note that the probabihty measures with finite support are weak 
star dense in A^(S"'~-^, 6). Since is continuous, we obtain, for all jj, G 



d-1 J 



C'^{{x,y))d^i{x)d^i{y)>Q, Vfc = 0,l,2,... 

1 

We can then estimate 

\{x,y)\Pdfiix)dfiiy) ^ [ [ f]A,C,"((x,2/))d^(x)dM(y) 



d-l ./ ed-1 



k = 

From the results in [2§|, one can deduce that 

l-3-5---(2t-l) 



Y^Xk [ I C^{{x,y))df,ix)d,xiy)>\o. 



An = 



+ 2) • • • (d + 2t - 2) ' 



which provides the desired estimate. 

We still have to address the "if and only if" part. Equality holds if and only 
if n satisfies 

C^{{x,y))d^lix)d^liy)=0, Vfc = l,...,p. 



d-l ./ Cd-1 



Wc shall follow the approach outlined in |34| in which the analog of Theorem 
was addressed for finite symmetric collections of points. In this case, the finite 
symmetric sets of points lead to finite sums rather than integrals as above. The 
key ideas that we need in order to use the approach presented in [s^l are: First, 
fl{E) := \{ii{E) + n{-E)), for E G B, satisfies PFP(/i,p) = PFP(/i,p). Thus, 
we can assume that fj. is symmetric. Secondly and more critically, the map 

\{x,yWdfi{x) 

is a polynomial in y. In fact, the integral resolves in the polynomial's coefficients. 
These two observations enable us to follow the lines in [34| , and we can conclude 
the proof. □ 
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Remark 4.11. One may speculate that Theorem 14.101 could be extended to 
p > 2 that arc not even integers. This is not true in general. For d = 2 
and p — 3, for instance, the equiangular FUNTF with 3 elements induces a 
smaller potential than the uniform distribution. The uniform distribution is a 
probabilistic tight 3-framc, but the equiangular FUNTF is not. 
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